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Abstract 

Starting from low energy effective chiral Lagrangian with gauged Wess-Zumino Wit- 
ten term, we have derived a hydrodynamic theory for chiral superfluid. It is a non- 
abelian hydrodynamics at zero temperature with only superfluid components. With 
an external electromagnetic field and baryonic and axial baryonic chemical potentials 
turned on, we are able to identify analogs of various anomaly induced term in normal 
hydrodynamics, including chiral vortical effect, chiral magnetic effect and chiral electric 
effect. As an example, we solved the hydrodynamic equations for the ground state and 
observed the chiral magnetic effect and chiral separation in the confined phase. 
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1 Introduction 



Successful applications of the fluid/gravity duality [1] to the R-charge black hole back- 
ground [2, 3] have found a vorticity induced contribution to the constitutive equation of 
the R-current. The vorticity induced term originating from a Chern-Simon term in the 
gravity description encodes chiral anomaly in the dual field theory. It signals that quantum 
effect like chiral anomaly can appear as parity odd terms in macroscopic hydrodynamics. 
It was later shown in a seminal paper by Son and Surowka [4] that parity odd terms are 
actually required by second law of thermodynamics for hydrodynamics with chiral anomaly: 
By demanding the existence of a entropy current with non-negative divergence, they were 
able to fix all the transport coefficients of the parity odd terms, which include chiral mag- 
netic effect (CME)[5] and chiral vortical effect (CVE)[6], in the constitutive equations of 
currents. Due to the non-renormalization of the anomaly coefficients, the anomaly induced 
effect are robust. There have been extensive model studies at both strong coupling using 
gauge/gravity duality [7, 8, 9, 10, 11, 12] and at weak coupling using field theory [13, 14]. 
Recently, the entropy current principle was also extended to higher dimension [15] and 
higher order hydrodynamics [16]. 

It was pointed out in [17] that the transport coefficients are fixed by the entropy 
current principle only up to unknown functions of temperature. It has been established 
later that the unknown functions are related to the gravitational anomaly [18]. The situa- 
tion becomes more complicated when the system contains spontaneous symmetry breaking. 
The Goldstone bosons resulting from the symmetry breaking is gapless, thus have to be 
introduced to the hydrodynamics as additional degrees of freedom, leading to superfluid 
hydrodynamics, the case of U(l) superfluid hydrodynamics has been studied in [19, 20, 21], 
based on early works on non-anomalous U(l) superfluid hydrodynamics [22, 23]. The dy- 
namics of the Goldstone boson enters the constitutive equations as superfluid velocity. 
The presence of the superfluid velocity lowers the predictive power of the entropy current 
principle, leaving several transport coefficients as arbitrary functions of thermodynamical 
quantities. 

There have been also effective field theory approaches [24, 25] to derive the consti- 
tutive equations of hydrodynamics. Since the anomaly arises as gauge variation of the 
Wess-Zumino-Witten (WZW) term, the effect of anomaly can be incorporated by introduc- 
ing a WZW term to the effective field theory action. One example is the chiral perturbation 
theory as the low energy effective field theory of QCD vacuum. The WZW term correspond- 
ing to the chiral anomaly in this case is explicitly known. The combined action can be used 
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to study the anomalous effect on the chiral superfluid hydrodynamics. This is essentially 
the approach taken in [24]. A different approach has been taken in [25], where the authors 
wrote down the effective field theory for fluid, and constructed the most general WZW 
term in dimension two. The resulting action allowed them to determine the anomalous 
transport coefficient explicitly. In this paper, we will follow the lines of [24]. We would like 
to emphasize here that since the chiral Lagrangian contains only Goldstone bosons as the 
fundamental field, what we will obtain is hydrodynamics with only superfluid components 
with non-abelian group structure. Our work distinguishes from [24] in many fine details 
as will be clear in what follows. These lead to quantitative differences in the resulting 
constitutive equations of the currents. 

This paper is organized as follows: In Section 2, we will focus on the abelian currents. 
We will derive the consistent current from the WZW term, and pass to covariant current 
by adding proper local terms in the external field. Section 3 will be a parallel study on the 
non-abelian currents. In Section 4, we will define properly the hydrodynamic variables and 
formulate a chiral superfluid hydrodynamics. Many interesting structures emerge from the 
WZW term, including the superfluid analogs of CME, CVE and also chiral electric effect 
(CEE)[21]. Section 5 includes an application of the hydrodynamics to the ground state of 
QCD, where we find the chiral magnetic effect and chiral separation effect. We will conclude 
in Section 6. 

2 Abelian currents from gauged WZW term 

Our starting point is chiral Lagrangian plus gauged WZW term[26, 27]. Since it is a low 
energy effective action, we will work in the mean field approximation, where all fields take 
their classical values. We will follow mostly the notations of [28] in this paper. The full 
action is given by: 



where U = e~ is the unitary matrix parameterized by the pseudoscalar <p and the pion 
decay constant /. (f> corresponds to the Goldstone boson originating from spontaneous 
chiral symmetry breaking. V = V^dx^ and A = A^dx^ are the one-forms of vector and 
axial vector flavor gauge fields. We further introduce the left and right one-forms from the 




(1) 



C = -LTr(D^UD»U + ), 



(2) 
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chiral field and the gauge fields: 

a = (d^U)U + dx^ = dUU + 
(3 = U + dU 

A L = l -{V + A) 

A R = \{V-A) 

In terms of these fields, the gauged WZW term is given by: 

T WZ (U,A L ,A R ) = T WZ (U) 

+ 5Ci J Tr(A L a 3 + A R (3 3 ) - 5C J Tr[(dA L A L + A L dA L )a + (dA R A R + A R dA R )f3] 

+ 5C J Tr[dA L dUA R U + - dA R d(U + )A L U] + 5C J Tr(A R U + A L U/3 2 - A L UA R U + a 2 ) 

+ ^ J Tr[(A L a) 2 - (A R (3) 2 ] + 5Cz J Tr{A\a + A%fi) 

+ 5Ci J Tr[(dA R A R + A R dA R )U + A L U - (dA L A L + A L dA L )UA R U + ] 

+ 5Ci f Tr(A L UA R U + A L a + A R U + A L UA R (3) 

+ 5C J Tr[A R U + A L U - A\UA R U + + ^(UA R U + A L ) 2 }. (3) 

The chiral field is SU(N) valued because in standard QCD, the diagonal U(1)a symmetry 
is broken by anomaly We will promote the chiral field to be U(N) valued as in [28]. This 
can be realized for example in the large N c limit, where the axial anomaly is suppressed. 
Under the general U(N)l x U(N) r gauge transformation parameterized by U(N) matrix 
valued ei and e R , 

5U = ie L U - iUe R 
5A L = de L + i[e L ,A L ] 
5A R = de R + i[e R ,A R ], 

C is left invariant, while the WZW term transforms as: 

5T WZ (U,A L ,A R ) = -WCi J Tr e L {{dA L ) 2 - \d{A\) \ - (L o R) . (4) 

The form of the WZW term in (3) treats the left and right chiral and gauge fields 
in a symmetric fashion, thus is not invariant under the either vector or axial vector gauge 
transformation. In reality, we need to impose invariance under the vector gauge transfor- 
mation, which requires subtraction of the Bardeen counter term [29]. The Bardeen counter 
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term corresponding to (3) has been worked out in [28]: 

r c = r wz (u = i,A L ,A R ). (5) 

Therefore, the complete WZW term is given by: 

T' WZ (U,A L ,A R ) = T WZ (U,A L ,A R )-T WZ (1,A L ,A R ). (6) 



The gauge transformation of the complete WZW term takes the following form: 



8T' WZ = 30Ci J Tr{e R - e L ) 



\Fv + + l^FvA 2 + AF V A + A 2 F V ) - l -A 4 



(7) 



with Fy = dV - {{V 2 + A 2 ) and F A = dA- \{V A + AV)* 1 . 

Now we wish to derive hydrodynamic currents from (1), (3) and (6). We require the 
hydrodynamic currents to be gauge invariant (covariant) for abelian (non-abelian) currents. 
In this section, we will focus on the abelian currents and generalize to the non-abelian cur- 
rents in the next section. For concreteness and possible connections to real world situation, 
we focus on three currents: baryonic current jg, axial baryonic current and electro- 
magnetic (EM) current jg. We will see that the baryonic current and the EM current are 
conserved while the axial baryonic current is not due to anomaly. We will turn on three 
external gauge fields: baryonic field B^, axial baryonic field B§ and electro-magnetic Aem, 
all in the Cartan subgroup of U{N)l x U(N) r . The first two do not exist in nature. We 
use them as a way to introduce baryonic and axial baryonic chemical potential. In the end, 
we put B^ = (n, 0) and B% = (/x 5 , 0). 

On general grounds, we expect energy-momentum (non)conservation to be one of the 
hydrodynamic equations: 



dpT-" = F^jl (8) 

where a = B,B§,Q corresponding to baryonic, axial baryonic and EM field respectively. 
The right hand side (rhs) is a sum over works done by external fields. We will see (8) 
naturally emerges as a result of the diffeomorphism invariance and gauge invariance. To 
derive (8), we start from the diffeomorphism invariance of the theory in curved spacetime: 

= SS = [ d 4 x4^Sg^ + [ d A x 5 4- T 8U + f d 4 x-^-SA a u . (9) 
J Sg^ J 5U J 6A« " w 

The diffeomorphism of the metric, chiral and external gauge fields are given by: 

S g ^ = -2V^V V \ 5U = V^d^U, 5A a = V u d u A a + {d^V u )A a u . (10) 



* lr There appears to be some minor sign errors in the original paper. We quote the corrected result 
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Note |jj = by equation of motion (EOM) of U, we are left with: 
= - J d 4 x2V»V u -^- + J d A x^(V v d u A^ + (dpV^Al) 

- j A^V^^) + J**V> (a„A^ - ( «L„-)) . (11) 

Since V u is arbitrary, we obtain: 

VT F = - Ald^\ (12) 

with the gauge invariant energy-momentum tensor = — —i= . Note that the current 
jn,a we obtained from the functional derivative is the consistent current. It is neither 
gauge invariant nor conserved. In order to formulate a hydrodynamic theory, we wish to 
have a gauge invariant current [17]. This can be done by expressing the non-conservation 
of J^ a as a correction to J^ a [25]: A u ' a d^ = -F£ a AJ»> a , so that the total current 

jH,a _ jfj,,a _|_ ^ jn,a - g g au g e invariant. 

To calculate (^jj^J, we consider the gauge transformation of the action: 

We can again drop the first term in the bracket by EOM of U . If the gauge symmetry 
is not spoiled by anomaly, SS = 0. For Cartan subgroup, we have in general 5A^ = d^e. 
Therefore we conclude that the current from functional derivative is conserved: 

If the gauge group is anomalous, 5S is nonvanishing, which should be used to obtain the 
nonvanishing We are interested in three particular gauge transformations list 

below, which correspond to EM gauge, baryonic gauge and axial baryonic gauge respectively: 

gauge Q: e L = e R = e, 5U = ie[Q, U], 5A em = 2de, 5B = 5B 5 = (15) 
gauge B: e L = e R = e, 5U = 0, 5B = 2de, SB 5 = 5A em = (16) 
gauge B5: e L = -e R = e, 5U = 2ieU, 5B 5 = 2de, SA em = 5B = 0. (17) 

They can be easily obtained from (4.17) of [28]. Note in the above, all e's are just numbers 
instead of matrices. 

From (7), we see Jb and Jq are obviously conserved, while Jb5 is not. We calculate 
its divergence as follows: 







SS 



J J B5 SB 5 = J 2J B5 de = -2 j ed,J B5 . (18) 
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On the other hand, our external gauge fields lead to 

V = B + A em Q, A = B 5 (19) 
F v = dB + dA em Q, F A = dB 5 . (20) 

Plugging them into (7) and compare with (18), we obtain 

15Cf ^Ci 
dJ B5 = — {(dBfTrl + 2dBdA em TrQ + (dA em ) 2 TrQ 2 ) + — (dB b ) 2 Trl. (21) 

Now we want to obtain explicit expressions for the consistent currents J R , Jb5 and 
Jq. We will focus on the contributions from the WZW term. The contributions from 
the original chiral Lagrangian can always be included later. It is easier to start with the 
non-abelian(matrix valued) left-right current as follows [30]: 

_ ST WZ (U,A L ,A R ) 
JL = 5A L 

= 5Cia 3 - 5C(2dA L a - A L da + 2adA L - duA L ) + 5CaA L a + 5Ci(-A L aA L ) 
+ 5C{-dUdA R U + + dUA R dU + + UdA R dU + ) + 5C{-U(3 2 A R U + - UA R U + a 2 ) 



+ 5Ci 



U{dA R A R + A R dA R )U + - 2dA L UA R U + + A L dUA R U + + A L UdA R U + 



- A L UA R dU + - 2UA R U + dA L + dUA R U + A L + UdA R U + A L - UA R dU + A L 
+ 5Ci(UA R U + A L a + aA L UA R U + - UA R {3A R U + ) 

+ 5C(A L UA R U + A L - UA R U + A L UA R U + ), (22) 
_ ST WZ (U,A L ,A R ) 

Jr = Mk 

= 5Ci(3 3 - 5C(2dA R p - A R d(3 + 2f3dA R - d(3A R ) + 5C(-/3A R /3) + 5Ci(-A R pA R ) 
+ 5C{-U + dA L dU + dU + dA L U - dU + A L dU) + 5C(U + A L Ufi 2 + U + a 2 A L U) 

+ 5Ci 2dA R U + A L U - A R dU + A L U - A R U + dA L U + A R U + A L dU + 2U + A L UdA R 

- dU + A L UA R - U + dA L UA R + U + A L dUA R + U + (dA L A L + A L dA L )U 
+ 5Ci(-U + A L aA L U + U + A L UA R /3 + (3A R U + A L U) 

+ 5C(-A R U + A L UA R + U + A L UA R U+A L U). (23) 
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The abelian currents we are interested in are simply given by: 

J B = \tt{Jl + Jr) - {U -> 1) (24) 
Jb 5 = \tt(Jl - Jr) - (U -> 1) (25) 
J Q = ^TrQ(J L + J R ))-(U-H). (26) 

Plugging A L = B + B ^+ A —Q and Ar = B-Bs+AemQ intQ ^22), (23) and (24), we obtain after 
tedious algebra: 

Jb = —Tr(a 3 + /3 3 ) - — (dBTr(a + p) + (L4 e m7VQ(a + /?) - A em TrQ(a 2 - /3 2 )) 
«/B5 = - — {dB 5 Tr(a + /3) + (L4 em TrQ(a - /?)) + — ^ em dA em (TrQ 2 - TrU+QUQ) 
J Q = — TrQ(a 3 + /3 3 ) + — 



- 3dBTrQ{a + 0) - dB 5 TrQ(a - /3) + B 5 TrQ(a 2 + /3 2 ) 



- (L4 em (2TrQ 2 (a + /3) + TrQdUQU + - TrQUQdU + ) + A em (TrQ 2 (a 2 - /3 2 ) - TrQdUQdU + ) 

— (dB 5 A em + 2dA em B 5 ){TrQUQU + - TrQ 2 ) - —B 5 A em Tr{QdUQU + + QdU + QU). 

(27) 

We can check explicitly that all three currents are invariant under B and Q gauges, but 
none is invariant under B§ gauge. The verification of the gauge invariance is lengthy. In 
the next section, we will see it is much more convenient to use the covariant variables. 

Now we want to correct (27) so that they are B§ gauge invariant also. We first write 
(12) out explicitly: 

= F%J* + + f*»j£ 6 - B 5X d,r m . (28) 

We wish to express B^dJsb as a sum of F * J terms. The former in terms of components 
of gauge fields, is given by: 



B 5\ d l* J B5 

T5Cz 

--D5A 



e^{d^B u d a B p Trl + 2d^B v 8 a A em ^TrQ + d^A em ^d a A em ^TrQ 2 ) 

(29) 



+ —e^d^B 5u d a B 5p Trl 



Indeed, it can be expressed as 

£ 5 ,A«5 = K AJ B + Fx^Jk + ^gAJg, (30) 
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where 

Aj£ = lhCie^d v B a B bp Tr\ + IbCie^ 'd v A em , a B bfS TrQ (31) 

Aj£ 5 = Wie^d v B ha B^Tr\ (32) 

A j£ = \^Cie^d v B a B^TrQ + IhCie^ d v A em ^ a B w TrQ 2 . (33) 

(30) can be proved by repeated use of the following identity: 

(d x X^ - d^X x )e^d u Y a Z p + (d x Y^ - d^Y x )e^d u X a Z p + Z x e^d^X u d a Y p = 0. 

(34) 

As a result, we obtain our final results of the gauge invariant currents: 

3b = Jb + AJ b , j B = Jb5 + AJ B5 , Jq = Jq + AJq. (35) 

We can check explicitly that the currents j R , 3b§ and jq are invariant under all three gauge 
transformation. 

3 Extension to non-abelian currents 

We also wish to obtain all non-abelian currents, which are (22) and (23) with properly 
chosen corrections to restore gauge covariance. We again turn on only the electro-magnetic 
field and the fictitious baryonic and axial baryonic fields. We do not turn on isospin chemical 
potentials, because in the chiral limit we are working at, a finite isospin chemical potential 
will lead to Bose-Einstein condensation of the associated Goldstone boson [31], thus further 
reducing the number of Goldstone bosons. 

It is convenient to work with the covariant quantities defined as follows: 

DU = dU - iA L U + ill A R , DU + = dU + - iA R U + + iU + A L 

a = DUU + , p = U + DU. (36) 

In terms of the covariant quantities, the consistent currents (22) and (23) take the following 
compact forms: 



J L = hCia 6 - 5C 



AidA L A L + 2(dA L a + adA L ) 



+ 5C(-DUdA R U + + UdA R DU* 



J R = 5Ci/3 3 - 5C 



- UdA R A R + 2{dA R p + pdA R ) 



+ 5C{-U + dA L DU + DU + dA L U). (37) 
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In arriving at (37), we have used A\ = A 2 R = and the commutative relation among Al, 
A R , dAi and dA R . Subtracting the Bardeen counter terms, we obtain: 

J c l = Jl-Jl{U^\) 

= 5Cia 3 - 10C(dA L a + adA L ) + 5C{-DUdA R U + + UdA R DU + ) 

- 5C(4idA L A L - AidA L A R - 2idA R A R + 2idA R A L ) (38) 

J c r = Jr- Jr(U -> l) 

= 5Cifi 3 - 10C(dA R p + j3dA R ) + 5C(-U + dA L DU + DU + dA L U) 

- 5C(—4idA R A R + 4idA R A L + 2idA L A L - 2idA L A R ). (39) 

It is straightforward to work out the gauge variations of the left and right currents. 
Starting with the gauge variations of the covariant quantities: 

gauge B: SU = 0, SA L = SA R = de, SDU = 5DU + = 0, da = Sfi = 0, (40) 
gauge B5: SU = 2ieU, 5A L = -5A R = de, SDU = 2ieDU, SDU + = -2ieDU + , 

5t\ = Sp = 0, (41) 
gauge Q: SU = ie[Q, U], SA L = SA R = deQ, SDU = ie[Q, DU], SDU + = ie[Q, DU + ], 

6a = ie[Q,a], 50 = ie[Qj], (42) 

we readily obtain the gauge variations of the consistent currents as follows: 

Q: SJt = ie[Q, J£], J C R = ie[Q, J C R ] (43) 
B: SJl = SJ c R = (44) 
B5: SJl = -5C(8idA L de + 4idA R de), SJ C R = -5C{4idA L de + 8idA R de). (45) 

Similar to the case of abelian currents in the previous section, the non-abelian currents 
are Q-gauge covariant and B-gauge invariant, but are not B5 gauge invariant. We need 
corrections to the currents to restore the gauge invariance. We note that the corrections 
are simply given by the U independent terms: 

AJ£ = 5C(4:idA L A L - AidA L A R + 2idA R A L - 2idA R A R ) 

AJ R = 5C(4idA R A L - 4idA R A R + 2idA L A L - 2idA L A R ). (46) 
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Our final results for the covariant currents are given by: 



f L = J c L + AJi 

= 5Cte 3 - WC(dA L a + adA L ) + 5C(-UA R U + A L UA R U + - DUdA R U + + UdA R DU + ) 

(47) 

f R = J c R + AJ R 

= 5C0 3 - 10C(dA R p + (3dA R ) + 5C{-U + dA L DU + DU + dA L U + U + A L UA R U + A L U). 

(48) 

We can verify that ^Tr(j£ +j R ), \Tr(j c L — j c R ) and ^TrQ(j£ + j c R ) reproduce the abelian 
currents j R , jB5 and jq in (35) respectively. 

The non (conservation) of the covariant currents jy and j c A can be obtained from (7). 
The most general gauge variation is given by: 



5S = j Tr 
= I Tr 



ss ss 

-{de L + i[e L , A L \) + —— (de R + i[e R , A R ]) 



I 
I 



5A L " ^ SA R 
e L (-dJ c L + i[A L , J£]) + e R {-dJ R + i[A R , J R ]) 

Tr(e L D L J c L + e R D R J R ), 



(49) 



where Dl = d — i[Al, ] and D R = d — i[A R , }. With our choice of the external gauge fields, 
the left-right covariant derivatives simplify to Dl = D R = d — ^f^lQ, ] = D®. Comparing 
with (7), we readily obtain 

D^Ji = W^Fy + 1^1), D®J R = -WCi{\F v + ^F\). (50) 

Combined with the corrections to the currents (46), we have 

D Q j c L = WCi(dA L ) 2 , D Q f R = -WCi(dA R ) 2 . (51) 



4 A hydrodynamic description of the chiral superfluid 
4.1 non-abelian hydrodynamics without anomaly 

In this section, we attempt to formulate a hydrodynamic description of the chiral superfluid. 
This is a zero temperature version of [32]. The new ingredients are the effect of anomaly and 
external electromagnetic field and chemical potentials. We start with the action without 
the WZW term. The stress tensor and currents can be easily obtained from functional 
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derivatives* 2 : 

e'2 



T ia , = --j {TrD^UD u U + - ^ u TrD x UD x U + ^j 

= ~ (tt (iD^UU + iD u UU + ) - X -g^Tr (iD x UU + iD x UU + ^ (52) 



where we have used an index in the currents to denote the contributions from the non- 
anomalous part of the action. The dynamical equations are simply 

d,T^ = FtfJ* + F%J* (55a) 
^ Vol = D Q ,4 R = 0, (55b) 

4q = h Tr Q {Jol + Jor) . <b 5 = \ Tr ( J ol ~ 4r) and F Q = dA™, F B5 = dB 5 . 
Note that the baryonic gauge field B is absent in (55a) and (55b). The axial gauge field 
i?5 is used to model the axial chemical potential B§ = (^5 , 0) . Semi-classically it originates 
from fluctuations of instantons in vacuum. We separate its contribution to the covariant 
derivative D^UU + = D^UU + - iB% due to its special role. 

In [32], the phase of the condensate U was chosen as the dynamical variables, we 
instead use iD® ^UU + and iU + JJ , which has an interpretation the non-abelian gener- 
alization of superfluid velocity. To justify this interpretation, we note that iD® ' yUU + and 
JJ are hermitian (In abelian theory, the hermiticity condition will be replaced by 
the reality condition for the superfluid velocity). Their decomposition on the basis of U(N) 
generators all have real entries. The hermiticity can be seen by noting the following identity 

D Q il UU + + UD Q fl U + = 0. (56) 

Therefore, we define the non-abelian superfluid velocity as: 

iD®^UU + = $ = ^ a T a , -iU+D^^U = tf = tf a T a , (57) 

with T a being the U(N) generators. In terms of the new variables, the constitutive equations 



# 2 The non-anomalous currents are automatically covariant and conserved. 
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become: 

JSl = t ^ + B ") 
4r = t - ^5 ) 

■% = yrrQ(£ L " + (58) 

Treating B^ as non-dynamical field, we can derive the conservation of energy-momentum 
tensor (55a) from the conservation of left and right currents (55b). Therefore (55a) is redun- 
dant thus can be discarded from the dynamical equations. We could make dynamical 
by writing B^ = fi^u^ and interpreting as a normal component in the hydrodynamic 
theory. Note that at zero temperature B^ originates from the fluctuation of instantons in 
vacuum. A dynamical B§ would mean we include the dynamics of the gluons. We will not 
do so in this paper, but only treat B£ as an external field. 

The defined superfluid velocities allow us to have a compact form for the constitutive 
equation. The price to pay is an increased number of unknowns. We have two 4— component 
superfluid velocities, in total 8 unknowns, but only 2 dynamical equations. The system of 
equations does not close. Note that components of & and £^ are not independent. It is 
not difficult to prove the following identities: 



D Q vi L _ jyQ^L = ^ + [Q> U]U+ 



D Q „$ - = ^] - U+[Q, U]^f. (59) 

We can close the system by taking the = 0,i component of (59), i.e. 

D Q t g _ D Q.^L = ^ _ {Q _ UQU+) El 

- DQ t tf = ef ] + (U+QU - Q)y, (60) 

which provide dynamical equations for ^ L and We have used = E{. Unlike the 
external magnetic field, which affects the dynamics, but not pump energy into the system, 
the external electric can possibly heat up the system and upset our zero temperature de- 
scription of the system. However, since the EM field remains entirely arbitrary, we can 
always have an electric field that is orthogonal to the currents and does no work to the 
system. 
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Note the appearance of UQU + and U + QU in (60) mean we have not complete the 
equations of motion. Defining the hermitian variables \ L = UQU + and \ R = U + QU, we 
easily find the equations satisfied by \ L an d X R '- 

D Q »x R = -i[$,X R \- (61) 

As in the case of £ LfJl and the time components of (61) provide the dynamical equations 
for \ L an d X R - To summarize, we have found the following hydrodynamic equations: 

D®[mJq L = D^^Jq R = 

D Q ti L _ D Q ti L = _ {Q _ ^ 

D^ R -D^ = i[^] + ( X ^-Q)^ 
D Q t X L = -i&,X L ] 

D Q tX R = -i[g,X R l (62) 

The first line provides the dynamical equations for £^ and The rest of the equations 
provide dynamical equations for \ L an d X R respectively. The additional equations 

given by (59) and (61) are non-dynamical. They are just constraint equations, which should 
be satisfied by the initial conditions. Once satisfied by the initial condition, the constraint 
equations will continue to hold as the system evolves. 

4.2 non-abelian hydrodynamics with anomaly 

Now we are ready to consider the effect of the WZW term on the hydrodynamics. The 
WZW term does not have correction to the LO of the stress tensor, which can be easily 
verified. It does induce corrections to the non-abelian currents. We have worked these out 
in the previous section. We quote the results as follows: 

j L = 5Cia 3 - WC{dA L a + adA L ) + 5C(-aUdA R U + - UdA R U + a) 

j R = 5Cip 3 - 10C{dA R p + j3dA R ) + hC{-U + dA L Up - pU + dA L U) (63) 
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In terms of explicit components, we have 



+ 



5Ci 



: fiuptr 



d u A e p m (2Q + X L ){£ + B 5a ) + + i?5 CT )(2Q + x L )d„A 



em 
P 



+ 



5Ci 



pupa 



d u A e p m (2Q + X R )(-^ + B 5a ) + (-£ + B 5a )(2Q + X R )d u A 



em 
P 



(64) 



where we have used the abelian nature of B£ to simplify the expressions. We can use (59) 
to rewrite (64) into a more suggestive form: 



5Ci 



+ 



5Ci 



.pvpa 



dvA™{^- + + B 5ct ) + (tf + + 



em 
P 



(65) 



- [^-(^ + ^)(c L - s 5ct ) + - B 5a )( 5 -Q + ^; 

Due to the effect of the WZW term, the hydrodynamic equations are slightly modified 
as follows: 

IhCi 



= MCi(dA L ) 2 = ^Pe^FQFQQ 2 
dQ 'Mr+ 3r) = -30Ci(dA R 



2 _ 15Cl a Q Q 2 
— o fc r pv r pa^ 



2 



(66) 



Note that the external baryonic field B^ or the baryon chemical potential does not appear 
in the hydrodynamic variables and equations. However, a nonvanishing baryon chemical 
potential can change the ground state of the theory, thus implicitly entering the expressions 
for the currents. We will see such an example in the next section. 

The baryonic and axial baryonic currents defined in (24) are of particular interest. 
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The relevant contributions are given by: 



Tr (4 L + 4 R ) = 

Tr {J£ L ~ 4 R ) = ^ (Tr(^ - + 2NBg) 



+ 5Cie pupa 



W v A™TrQ(t L a - g) + 6d u A e ™TrQB 5a 



Tr (f L ~ f R ) = 10Cie!«"»d v A™TrQ(g + g). 
The anomalous baryonic current can be written alternatively as: 
Tr (f L + f R ) = -lOCe^Trgtfg 



(67) 



+ 5Cie^ pa 



3d v A™TrQ(g - g) + Qd v A™TrQB ba 



(68) 



which agrees with the previous results [30, 33] upon setting B p = 0. It is illuminating to 
separate the temporal and spatial components of the currents induced by WZW term: 

Tr(f L + f R ) = -hCi^ k Tr (D^^k ~ D Q ^) + l5CiB t TrQ(tf - g) 
Tr(f L +j R ) = -5CiTr(nf-nf) 



+ IhCi 



-R 



Tr(j° L - j° R ) = 10CiB t TrQ(tf + g) 

Tr(f L - f R ) = 10Ct (j^EjTrQitf + " B t TrQ(^ + & 



(69) 



where we have used B p = Qit 5 ,0). We have also defined nf = e ivpa TrD v i L p £ and Vtf = 
e tupcr TrD u ^^, which are analogs of vorticity in non-abelian superfluid, thus giving rise 
to the CVE. They are vanishing for abelian superfluid. This can be easily shown by using 
(59) and the commutativity of abelian superfluid velocities. Ei and Bi in (69) are external 
electric and magnetic field. The latter is not to be confused with components of B§, which 
originates from the fluctuation of instantons. The terms proportional to Ei and Bi can be 
interpreted as the analogs of the CEE and CME. Note that CVE and CME have already 
been found in [24]. The identification of CEE is new in the current analysis. From various 
terms in (69), we observe that there is CVE in the baryonic current and CEE and CME 
are present in both baryonic and axial baryonic currents. 
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5 Application to the ground state 



As an example, we calculate the currents in the ground state of QCD at finite baryonic or 
axial chemical potential in the presence of a constant external magnetic field. This will be 
done by solving the hydrodynamic equations (66) for the superfluid velocities, the results 
of which will be used to obtain the baryonic and axial currents from (67). 

We focus on one-flavor case in which we have the diagonal r/ meson only. The 
calculation simplifies significantly as all the non-abelian structures drop out. It is convenient 
to use the a and ft variables in solving the hydrodynamic equations (66). Since the field U 
is just a U(l) phase: U = e t,p , we obviously have a = ft, x L = X R = Qi therefore \ L an d 
X R do not enter the hydrodynamic equations. (66) reduces to the following: 

W8 + $) = 0, (70) 
where 

Jl = z^j-a , Jr = 

j L = -15CdA em a ,j R = -15CdA em a. (71) 

We proceed further by considering a plane wave ansatz: U = e %k>iX>i , thus we have = 
i(k^ — -B5 ), with B§ = (//5, 0, 0, 0). The constant magnetic field is applied along the z-axis: 
B = (0,0, B). It is easy to see all the equations are trivially satisfied. To determine the 
momentum k^, we need to minimize the free energy H — /xNb and H — fi^Nsb f° r finite \x 
and ^5 respectively [34]. Let us consider the case \i / 0, ^5 = first: 

H - fiN B = J d 3 x (J^- (jfcg + k\ + kl + kl) - IhCiuBk^j . (72) 

It is obviously minimized at k^ = (0, 0, 0, — 60 ). The resulting nonzero components of 
the baryonic and axial currents are given by: 



jo = -v ; ^ = 15CifiB . (73) 



{30CrffiB 2 
T 



Similarly for the case \x = 0, ^5 7^ 0, the free energy is given by 

H - » 5 N B5 = J d 3 x (j((ko + W;) 2 + kj + k 2 2 + kl) - ^ 5 (k + Ms)) • (74) 

It is minimized at k^ = (0,0,0,0), with the following nonzero components of the baryonic 
and axial currents: 

4 = 15Cifi 5 B, J° m = Pi!*. (75) 
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We identify in (73) and (75) the CME and chiral separation [35, 36] in confined phase. 
Also a finite baryonic and axial chemical potential induces corresponding charge density 
respectively, although they differ in the dependences on the external magnetic field. 

6 Conclusion 

Following the idea of [24], we have derived, at the mean field level, a non-abelian superfluid 
hydrodynamics from the chiral Lagrangian and WZW term. It is a zero temperature hydro- 
dynamic theory, with the dynamical equations involving the non-abelian currents only. The 
constitutive equations for the currents contain contributions from the chiral Lagrangian and 
the WZW term. We identify superfluid analogs of CVE, CME and CEE in the contributions 
from the WZW term. The hydrodynamic equations can be useful in the study of the low 
energy dynamics of QCD. As an example, we solved the ground state of the hydrodynamic 
theory at finite baryonic or axial chemical potential in the presence of a constant magnetic 
field. We obtained the anomaly induced contributions to the baryonic and axial currents, 
which can be viewed as CME and chiral separation in the confined phase. 

Throughout the paper, we have been restricting ourselves to topologically trivial con- 
figuration of the chiral field. It has been known that topologically non-trivial configuration 
such as vortex and domain wall allows for more terms in the effective Lagrangian [35, 33]. It 
is interesting to investigate the anomaly induced effect with a relaxed topological condition. 

It is also interesting to consider the effect of mass term and isospin chemical potential 
for the chiral field, which can lead to the Bose-Einstein condensation of the associated 
Goldstone boson. It is curious to see how the interplay of chiral condensate and BEC will 
modify the hydrodynamic description. 
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